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We describe the stress energy of a fluid with two unequal stresses and heat flow 

in terms of two perfect fluid components. The description is in terms of the fluid 

■ velocity overlap of the components, and makes no assumptions about the equations of 

^ . state of the perfect fluids. The description is applied to the metrics of a conformally 

O , 

, flat system and a black string. 

in 

PACS numbers: 95.30.Sf, 04.60.Sf 

I. INTRODUCTION 

Descriptions of fluid systems used in general relativity have become more complex as we 
acquire new data about our universe. Simple perfect fluids are still useful models in many 
cases p|, but, increasingly, extensions of the perfect fluid stress-energy to include anisotropy 
and fluid interactions are necessary considerations. The stress-energy for an anisotropic fluid 
can easily be written in terms of metric based tetrads [Ua, Ra,Qa, ^a] 

Tab = {e)UaUb + {Pr)RaRb + {Pe)<daeb + {P<p)K^b. (l) 

o3 ' without addressing the physical origins of the anisotropy. Fluid anisotropy has many physical 
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For a low density fluid, anisotropy can be modeled by multi-perfect fluid 



a 



or equations of state 



descriptions with differences in the component fluid velocities 
generating the anisotropy. One method of incorporating these differences is through stress- 
energy equivalence [g], Q] where a composite system can be written as the sum of the two 
component perfect fluids: 

{e)UaUbHPr)RaRbHPe)<da&bHP^)^a^b = {ei+PlM^U^^^ 

(2) 

This form allows the component fluid motions and individual equations of state to be built 
into the composite fluid description and has parameter freedom for modeling assumptions. 



A related method is due to Letelier [8|, where he exphcitly rotates the component fluid 
velocities to create a new composite tetrad. Using the new tetrad, the composite stress- 
energy tensor describes a fluid with one anisotropic stress, a, associated with a direction of 
fluid anisotropy T^. 



Tab = {e + U)UaUb + {cr- n)t„t, + {U)gab 



(3) 



Letelier's stress-energy form follows from an assumption of zero heat flow in the composite 
fluid and implies a relation between the stress-energies of the component fluids, ei + Pi = 
const X (£2 + -P2)- Not all physical examples will obey this stress-energy relation or will have 
zero heat flow, and Eq.(l3]) can follow from other physical assumptions. Multi-fluid models 
of complex fluids are increasingly being used, not only in formal general relativity [9] but 
also in other physics sub areas. Two and three fluid descriptions cover phenomena like multi 
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cosmological models 



171], particles in heavy ion collisions [18|, Fermi-Bose fluids [19| and hydrodynamic nuclear 
models 20j. A need has emerged for a range of two fluid descriptions to use as a basis 
for n-fluid generalizations 2ll425l]. In this paper we suggest an anisotropic two fluid stress 
form based on the overlap of the component fluid velocity vectors rather than a stress- 
energy assumption. The description allows heat flow along the direction of anisotropy in the 
composite stress-energy and has no initial stress-energy assumptions. In the next section 
we set up the stress-energy description. The relation between the 4-velocity overlap and the 
stress-energy is discussed in the third section, and some examples are discussed in the fourth 
part of the paper. 



II. COMBINING COMPONENT FLUIDS 



The stress-energy form 

Consider a manifold which contains two perfect fluid flows and metric gab- A composite 
stress-energy for two perfect fluids is 

Tab = (ei + PM'^UI'^ + {62 + P2)U^M^^ + gab{Pl + ^2)- (4) 

In order to express the composite stress-energy in a simpler form, following Letelier jsl, we 
map the timelike unit vectors \U'^\ U^^] to an un-normed pair [?7*, T*] where U* is timelike 
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and T* is spacelike. A general unimodular transformation between these two sets of vectors 
can be written as 

A cos a —B sin a 
D sin a C cos a 

with unit transformation determinant AC cos^ a + BD sin^ a = 1 and transformation inverse 









a 



(5) 



(6) 



C cos a B sin a 
— Dsina A cos a 

Since U*"' and T*" are not normalized, two of the constants can be absorbed into their 
definitions. Choosing A = C = 1, the unit determinant condition is 



cos a + BD sin a = 1 



(7) 



The rotation angle a is fixed by requiring g"'^T*JJl = 0, with U* timelike, T* spacelike and 

-1. The orthogonality condition implies 



{B - D) cos a sin a + f/^^^'^f/i^) (cos^ a-BD sin^ a) = 0, 



(8) 



where f/'-^-'^f/i^'' is the velocity overlap. There are several special values of B which can be 
eliminated from the parameter range by the requirement that T* be spacelike. Condition 
= in determinant Eq.Q implies a = 0, vr. From the transformation equation Eq.([6]), this 
choice also identifies both U* and T* as timelike vectors so that overlap equation ([8]) is not 
valid for values B = 0, a = 0,tt. The choice B = ±1 requires D = ±1. The overlap equation 
thus becomes cos2q; = 0, with a = tt/A, Stc/A. These values of a in cos2q; = are 

also excluded by the spacelike condition on T*. From Eq.(|6]) with i? = ±1, we have 



T* = T sin aUi^^ + cos aU^^^ 
TIT*"" = -1 ±sin2«[-f/(i2)] 



(9) 
(10) 

(11) 



The velocity overlap U'^^''^^ must be negative, since it is the product of two future pointing 
timelike vectors. For T* to be spacelike requires ± sin 

2a[-t/(i2)] > 1, or a< 7r/4 for 5 > 0. 
The ranges we will consider are -B 7^ 0, 1 and < a < tc /A. These ranges and the unit 
determinant condition require BD = 1. The two normalized unit vectors are 



U" 



(12) 
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with norms 
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(13) 



III. OVERLAP AND STRESS-ENERGY 



Stress-Energy 



The rotation angle can be expressed in terms of the 4-velocity overlap. 

, 2B 



tan 2a 



1,2) 



B^0,1 < a < 7r/4 
The composite stress-energy tensor in Eq.(j4]) can now be written as 

Tab = {e + U)UaUb + {cr- n)t„t, + {U)gab + QaUb + QbUa 
with fluid parameters 

U=pi+p2 

Qa = f aN N [B~\e2 + P2) - B{e I + Pi)] sin a COS a 
€ + U = N'^[{£i + Pi) cos^ a + B-^{e2 + P2) sin^ a] 
a - n = N'^[B'^{ei + Pi) sin^ a + (£2 + P2) cos^ a] 
and B an unspecifled constant except for P 7^ 0, 1. 



(14) 



(15) 

(16a) 
(16b) 
(16c) 
(16d) 



Heat Flow and B 

Letelier's {s] choice, P^ = (£2 +P2)/(^i + Pi) reproduces the original two fluid stress- 
energy with Qa = 0. However, the underlying physics of the description can depend on both 
the velocity overlap and general equations of state in the component fluids. For example, 
if the two fluids move together with the same 4-velocity, the overlap is —1 and we have 
tan(2ao) = 2P/(1 — p2). Another way of choosing P that allows non zero heat flow is to 
use the f/i^^ = Ua^^ condition to deflne P as 



P = tan ao 



(17) 
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with a = ao producing the 'ahgned' fluid. Substituting the ahgnment condition in Eq.(|9]) 
requires the ahgned fluid to have zero anisotropy vector, T* = [recaU, U* = cos aUa^^ + 
B sin aUa^^ , T* = (— sin aUa^^ + cos aUa^^] . The conditions on B set the positive range 
< ao < it/A, such that ao < a < n / 4. The velocity overlap is 



tan 2a = -^/^^'^^ tan2ao 



(18) 



The general normalizations are 



1 



2 sin ao 
1 



cos 2ao 
cos 2a 
cos 2ao 
cos 2a 



1 

+ 1 



(19) 
(20) 



2 cos^ ao 

Note that both and T* have zero values when the fluid is aligned with a = ao. The fluid 
parameters are 



n = Pi + P2 



sin 2ao 



sin 2ao 



cos2ao 
cos 2a 
cos2ao 
cos 2a 



+ 1 
- 1 



[(^i + Pi) sin^ ao cos^ a + (£2 + P2) cos^ ao sin^ a] 
[(^i + Pi) sin ao sin a + (£2 + P2) cos ao cos a\ 



(21a) 
(21b) 

(21c) 



The heat flow vector is 

/ sin a cos a / cos 2ao 
\ 2sin2ao V^^^^ 

and is zero for a = ao. 



|2 - 1 [{62 + P2) - tan^ ao(£i + Pi)] 



(22) 



The a, ao boundary 

The tf'^a^ = Ua^^ alignment condition has stress-energy relations Q = and e + H = 
El + Pi + 62 + P2- When a is close to ao, the 4- velocities should be only slightly different 
and the fluids should show small deviations from the aligned relations. There are three 
parameter boundary cases: (1) both a and ao near zero, (2) both a and ao near n/A, and 
(3) ao near zero and a just under 7r/4. 

One expects the flrst two cases to describe a composite fluid only slightly different from 
the ahgned composite. For the flrst case, ao = So, a = 6, choose ao and a small but of the 
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same order with Sq < 5. We have R ^ 1 + 25'^ -25^ and tan2a/VR'^ - 1 ~ 5/ ^J5'^ - 51 » 1. 
The fluid parameters are 

-f/(i2) _ ^/^^ ^23a) 

£ + n ~ (£i + Pi) + (^2 + P2)(<5/<5o)' (23b) 
a-n~ [(5/5o)'-l](^2 + P2) (23c) 

with heat flow 

Q ~ (<5/5o)v/(V5o)=^-l (^2 + P2) 
The second case, a = n/A — 6, ao = it /A — 60, is very similar to the first with 60 > S. The 
velocity overlap is again, almost aligned, — f/^^'^^ ~ ^o/^^ but the composite fluid relations 
to the component fluids are multiples of the aligned fluid description. 

g ^ 1 - 1 [{62 + P2 - (£1 + Pi)] (24a) 

£ + n ~ ^ [60/6 + 1] [si + Pi + £2 + ^"2] (24b) 

a-n~ i[(5o/'5-l][£i + Pi+£2 + P2] (24c) 

The stress-energy in both these cases obeys the Ua^^ = Ua^^ condition to lowest order, 
e — a + 2n ~ £1 + Pi + £2 + -^2- The physical difference between the two cases can be 
explained by the composite * vector relation to the component fluid velocities, Eq.(j6]). For 
parameter values near zero, the component fluid velocity is dominated by the first fluid, 
U* ~ ui^^ with direction of anisotropy T* ~ —{6/6o)u!i^^ + Ua^\ For parameter values near 
7r/4, U* ~ {U^^^ + UP)/V2 and ~ {-U^ + uP)/V2. The straight combination of 
velocities in the second case explaining its close similarity to the aligned fluid case. The 
third case is an example of strong non-alignment. When a and ao are at opposite ends 
of the parameter range, ao = <^o, « = tc/A — 6, we have R = cos25o/ sin 25 >> 1 and 
tan(7r/2 — 25) j^R? — \ ~ 1. The fluid parameters obey a very different relation than the 
aligned fluid with e -|- 211 — cr ~ 0, rather than the sum of the component stress-energy. We 
also have 



(^2 +jP2 
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For this case the * vectors are related to the component velocities by U* ~ {Ua^^ + 6oUa^^)/ V2 
and ~ (-t/i^V'^o + Ua^^)/V2. In the next section we give some metric examples. 

IV. APPLICATIONS 

As an application of the two perfect fluid description we consider three different examples. 
The first two are metric based with an anisotropic stress-energy following from the field 
equations. The inverse of Eqs.fl21bJ?Ib) give the component stress-energies in terms of the 
composite descriptions: 



El + Pi = cos^ ao 



62 + P2 = sin^ ao 



cos 2a + 1 , cos 2a — 1 

ie + n + (T - n 

cos 2ao + cos 2a cos 2ao — cos 2a 

cos 2a — 1 , ^, cos 2a + 1 
{e + n) — + {a- n)- 



(26) 
(27) 



cos 2ao + cos 2a cos 2ao — cos 2a 

The two metric examples are a conformally fiat spacetime and a black string. The stress- 
energy from the field equations for both spacetimes has e + 11 = 0, and we can describe the 
related component perfect fluids. The third example uses two dusty component perfect fluids 
and a single anisotropic stress with no heat flow. The component fluid density relations are 
examined along with the equation of state in the composite. For this example, the Letelier 
description and the description in this paper coincide. 

Example: A conformally flat spacetime 

A simple conformally flat spacetime has metric and field generated fluid parameters as 
seen by a comoving observer f/" = (e~"^, 0, 0, 0). 

ds"^ = e^'^'i-dt^ + dr^ + r^d(j)^ + dz^), (28a) 
8716 = -a2e-2'^^ (28b) 



in 



n = a'e-'"', (28c) 



U = Pr= P^, (28d) 

8na = SirP^ = ^a^e-^""', (28e) 

with the anisotropy in the z-direction, = (0, 0, 0, e~°^). The z-dependent negative density 
in this solution does not lend itself to a physical description outside of a cosmological constant 
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or Casimir effects. The two ffuid description of the composite ffuid model has the advantage 
of explaining the negative composite density in terms of tension in the first fluid. For this 
fluid £ + n = and 87r(cr — 11) = 2a^e~^'*^. From Eq. (l2T]) . the component fluid parameters 
that could create a composite fluid are 

2a^e~^''^ (cos 2ol — 1) cos^ 

^1 + "i = — 5 7 7. (29) 

87r (cos zao — cos la) 

20^6"^""^ (cos 2« + 1) sin^ 

£2 + i^2 = — 7 ^ 30 

87r (cos 2ao — cos la) 

The second fluid could have both positive stress and density. The first fluid, if it has positive 
density, must describe a fluid with tension rather than pressure. The heat flow is axial 



(c?e 2a2gjj^2a /cos2ao + cos 2a \ ^ , , 

\ 87r V cos 2ao ~ cos la j 

The 4- velocities of the component fluid come directly from the transformation equations and 
involve, unsurprisingly, a coordinate velocity in the z-direction. 

f/(^)" = [e-'^^A^cosa, 0, 0, - e"""iVtanaosina] (32) 

^ rg-a.^i^^ C-'^^iV COS d (33) 

tanao 

was imposed in Eq. (ll2p . Both U^^^°- and f/^^)'' are expanding, accelerating 
and shear-free. Requiring the velocity overlap, Eq.f llOp to be constant, f/*^^^) = const, can 
be restrictive for shear- free fluids, requiring acceleration and expansion (see Appendix A). 
That restriction is met for this example. 



Example: A black string 



A similar example to the conformal metric is the black string of Lemos and Zanchin 26| . 
with metric 

ds^ = -(Ar^ - m)df + + r2rf02 + rf^2 r^^^ 

yhr^ — m) 

and comoving stress-energy following from the field equations 

Svre = -A, (35a) 
Svrn = A, (35b) 
Sttct = SvrP^ = 3A. (35c) 
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From the first example, the replacement a^e~^"^ — )■ A, gives the component stress-energy and 
heat flow for this case. The direction of anisotropy is along the string axis. Here, since the 
cosmological constant can be negative and does not depend on position, a negative energy 
density is possible. The fluid matter obeys e + 11 = and a = 311 . The 4-velocities for the 
black string components are 

N 



jj{2)a 



A/Ar^ — m 



N 



cos a, 0, 0, — iV tan ao sin a 



sm a ~ , 

I , , 0, 0, A^cosd 

VAr2 - m tan ao 



(36) 
(37) 



These component velocities are expansion-free but have acceleration, shear, and vorticity. 

Example: Linear composite equation of state 

For the third example consider a composite fluid with 11 = 0, and two component dust 
fluids. Pi = P2 = 0. The composite fluid has a linear equation of state related to a and ao'. 



a 



cos 2ao I 
cos 2a 



cos 2a() 

cos 2a 



(38) 



With Eq. (ET|) we have for the composite fluid 



El sin^ ao = 62 cos^ 



(39) 



and this example has the Letelier [8| stress-energy relation with zero heat flow. 



£1 



a 



2 cos^ ao 
2 cos^ ao 



cos 2ao 
cos2q; 
cos 2ao 
cos2q; 



(40) 
(41) 



The composite fluid has a dusty equation of state for a close to ao, R = (cos 2ao/ cos 2a) ~ 1, 
changing to a stiff equation of state for large R. 



V. DISCUSSION 



In conclusion, we have presented an anisotropic fluid recipe based on component velocity 
overlap. It is useful in describing anisotropic fluids in terms of possible perfect fluid compo- 
nents, and may also be of use in developing broader descriptions of multi-fluid models. The 
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description includes heat flow driven by fluid velocity non-alignment. The heat flow is along 
the direction of anisotropy in the stress-energy form. Heat flows with several spatial vector 
components are possible [9], especially if spatial tetrads are chosen to describe acceleration, 
or a null-vector construction, or other non-heat related parameters. Including more general 
heat descriptions is an idea for future work. The fluids considered here are non-interacting; 
another possible extension is to component fluids with non zero, but balancing stress-energy 



divergence 



32-[35|], and to acoustic phenomena related to the heat flow |36l-l38|. 



Appendix A: Restrictions on the Component Fluids 

The behavior of the composite fluid is determined by the component velocity overlap and 
equation of state. The velocity overlap equation following from tetrad orthogonality can 
be restrictive since it requires the velocity overlap to be constant, f/^^'^) = const. This 
will impose restrictions on the component fluids. Consider the covariant derivative of the 
constant velocity overlap Vfcf/^^'^-* = 

tjiDa^^) = -u^^'^um (Ai) 

Each of the velocity derivatives can be expanded in terms of its acceleration, expansion, 
shear and vorticity {k\ 6*, cr*^, w*^) and projection operator h^^ = gab + Ua^'uf^ . 

The fluid acceleration can be parameterized with the Frenet tetrad associated with each 
velocity vector [Ua \ Aa\ Ba \ Ca^]- The acceleration lies along the vector Aa^ and can be 
written as Ifa^ = U^\u^''^^ = K,Aa\ The acceleration for fluid one can be isolated with t/^^^'' 
multiplication and similarly, for fluid two with U^'^^'^ multiplication: 

3 

^^^(2)f>(l)a ^ ^^i(l)f/(2)af;{l,2) _ ^^^^f/{2)af>(2)6 _ + (f/(1.2))2] 

3 

If the component fluids are shear-free we have 

3 

^^^(2)f>(l)a ^ ^^i{l)f/{2)af;(l,2) _ ^ (f/(l'2))2] (A4) 

o 
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If the two component fluids are shear free and have no expansion then they satisfy 

If the expansions are zero, this can only be satisfied by identical component fluid velocities, 
f/^^'^) = —1 or by zero accelerations, ki, K2 = 0. The component fluids must be unaccelerated 
or aligned. If the fluids are not aligned or if there is acceleration, one (or both) must 
have expansion if they are shear-free. In the first metric example, the component velocities 
Eqs. fl32|33l) are shear and vorticity free, and have both expansion and acceleration. The 
black string component fluids have zero expansion but non-zero radial acceleration, with 

ri 

shear components aor and azr, and non-zero vorticity. Collins [27] has conjectured that, 
if shear-free perfect fluids obey a barotopic equation of state, then either the expansion 



or vorticity must be zero. If the conjecture is true |28l-l31l|. the shear- free, accelerated, 
unaligned component fluids must have zero vorticity or a non-barotopic equation of state. 
The conjecture holds for the conformal metric and does not apply to the black string. 
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